NACA TN 3069 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 3069 




INCOMPRESSIBLE FLOW PAST A SINUSOIDAL 
WALL OF FINITE AMPLITUDE 
By Carl Kaplan 

Langley Aeronautical Laboratory 
Langley Field, Va. 



MC 

£Fi 251! 


Washington 
February 1954 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 


NATIONAL ADVISORY COMMITTEE FCR AERONAUTICS 


TECHNICAL NOTE 3069 



tsia 


IHCOMPRESSIBIE FLOW PAST A SINUSOIDAL 
WALL OF FINITE AMPLITUDE 
By Carl Kaplan 

SUMMARY 


The problem of the present paper has been chosen for the purpose 
of exhibiting some of the possible mathematical troubles that may occur 
in the iterative procedures so much in use in present-day aerodynamics. 
First, it is shown that the example of incompressible flow past a sinus- 
oidal wall of finite amplitude should properly be treated in the plane 
of the velocity potential and the stream function rather than in the 
physical flow plane. Then, two contrasting iterative procedures are 
utilized for the solution of this particular problem. One is the well- 
known small -disturbance method in which the physical-plane coordinates 
are determined in the form of Fourier series whose coefficients are 
a n alytically developed as series in ascending powers of the amplitude. 

In general, this method precludes ary discussion of convergence, the 
tacit assumption being that no mathematical limitation intervenes before 
the solution ceases to be valid because of some physical reason. The 
other procedure is to state the problem in the form of an integral equa- 
tion whose solution can be found by a process of successive approxima- 
tions. The convergence of this method can usually be judged when the 
difference between any two successive approximations is deemed negli- 
gible. An included numerical example serves to emphasize the superior- 
ity of the integral-equation approach over the small-disturbance method. 


INTRODUCTION 


In the nonlinear treatment of stationary compressible flows, the 
methods utilized very often involve developments in series and succes- 
sive approximations. One of the most widely used of these methods is 
the Prandtl-Busemann small-disturbance procedure for the calc ulat ion of 
two-dimensional compressible flow past a thin profile . According to this 
method, the velocity potential or the stream function is developed in 
powers of the thickness coefficient with the thin profile placed at a 
vanishingly small angle of attack in a uniform stream of Mach n umb er less 
than unity. Thus far, no strictly mathematical investigation has been 
made on the convergence of this small-disturbance procedure. Indeed, 
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•the labor increases so rapidly with the order of the approx im ation that 
usually only two or three steps can be calculated; thus, any possible 
rigorous statement on the question of convergence is eliminated. The 
first step in this method of iteration being the undisturbed stream 
itself, the decisive element in the convergence of the process is the 
tacit assumption that no purely mathematical limitation sets in before 
the breakdown of the flow due to some physical reason (for example, the 
attainment of sonic velocity at some point in the field of flow) . That 
such mathematical difficulties can appear, even without the complication 
of compressibility, is illustrated in the problem treated in the present 
paper, namely, incompressible flow past a sinusoidal wall of finite 
amplitude . This solid boundary has been chosen chiefly because the basic 
assumptions of the small-disturbance method are adhered to, there being 
no stagnation points or points of infinite velocity in the field of flow, 
and because numerous iteration steps can be achieved without undue labor. 
Moreover, it was found possible to solve the problem independently by 
means of an integral-equation approach much in the manner of the arbitrary- 
airfoil theory of Theodorsen and Garrick (ref. l). A comparison of these 
two modes of solution clearly reveals basic weaknesses in the small- 
disturbance method when it is utilized for the purpose of approximating 
analytically the exact solution of a flow problem. 

It may be mentioned that the sinusoidal wall in two dimensions and 
its axisymmetric counterpart, the corrugated cylinder, have in the past 
been of considerable aid in the solution of various problems in aero- 
dynamics (refs. 2 to 7 ). It is believed that the material of the present 
paper may be found useful in the treatment of such diverse problems as 
panel flutter and the hydrodynamic theory of water waves of finite 
amplitude . 


AMLISIS 
General Formulas 


When the Prandt 1-Bus emann small-disturbance method is utilized to 
obtain the two-dimensional incompressible flow past the wavy wall 

y = e cos x 

the complex potential of the flow w can be represented in the following 
fashion: 



Bn(e)e 


inz 


( 1 ) 
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where 

x,y nondimensional rectangular Cartesian coordinates, radians 

z = x + iy 

e amplitude of wavy wall, radians 

0 nondimensional velocity potential 

\Jr nondimensional stream function 

B n (e) real quantities presumed to be expansible in ascending powers 
of e 

and where U, the undisturbed stream velocity, and A, the wave length 
in radians, are utilized as units of velocity and length, respectively. 

According to equation (l), it is tacitly assumed that w is an 
analytic function of z along and above the solid wavy wall. However, 
the covergence of the series appearing on the right-hand side is con- 
fined to the broadest strip which is parallel to the x-axis and is 
devoid of singular points of w (ref. 8). Thus, in order for this 
series to represent w on and above the sinusoidal boundary, w must 
be free of singularities in a strip containing the boundary. There is 
no assurance, however, that such singularities do not appear in the 
portion of the strip below the wavy wall as soon as the solid boundary 
departs from a straight line. Therefore, the solution in the form of 
equation (l) is not always valid and may even be divergent for e > 0. 

The rather surprising statement that the small-disturbance itera- 
tion method cannot always be applied in the physical plane to the pro- 
blem of flow past a sinusoidal boundary poses the question of what is 
the correct mathematical approach to this well-defined flow problem. 

The answer is simply that the problem should be treated in the 0ijr-plane, 
where the sinusoidal boundary y = e cos x is replaced by the straight 
line = 0. Equation (l) is then replaced by 


z = w + i iAo(e) + i f An(e)e ilW (2) 

d n=l 

where the right-hand side represents a complex function of w which is 
now analytic in the entire upper half -plane . The coefficients An(e) 
are reel quantities presumed to be expansible in ascending powers of e. 
No boundary condition need be stated in the ±0 directions because the 
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solution is periodic in 0; whereas for \|r = °°, the complex velocity is, 

as it must he, ^ = 1. Expressed in real form, equation (2) is equiva- 
dz 

lent to the pair of equations 

x = 0 - A Q (e)e“ ni l r sin n 0 (3a) 

n=l 


and 


y = if + i A q(s) + y Aix(e)e“ n ^cos n0 


n=J 


At the solid boundary, i}r = 0; therefore. 


(5b) 


x = 0 - Aj^eJsin n0 
n=I 


(4a) 


and 


e cos x 


0P n 

= \ Ao(e) + 2_ An(e)cos n 0 


n=l 


(4b) 


From equation (4b) (since e cos x is an even periodic function of 0), 


An(e) 


nTt 

-2 e f 

* J 0 


cos x cos n0 d0 (n = 0, 1, 2, . . . ») 


Inserting the expression for x given by equation (4a) into this equa- 
tion yields the following formula: 


*“ (£ > “f e X cos n0 cos 


°2_ 


0 - 2_ Am(e)sin m0 


m=. 


d0 


(5) 
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If the quantities A n (e) axe assumed to be expansible in a series of 
ascending powers of e, equation ( 5 ) represents the recursion formula 
for the determination of the coefficients of these series. Preliminary 
calc u l a tions show that the An(e) can be expressed as 


\ AqC 6 ) ~ y" a Om e2m (6a) 


and 


An( e ) 



a m e n+2m 


(n = 1, 2, . . . co) (6b) 


where the a ^m are numerical coefficients uniquely determined by 
repeated use of equation ( 5 ). Thus, the initial step in the determina- 
tion of the a nm is to obtain a-^Q, the coefficient of the first term 
of the series for Aq(e). From equation ( 5 ), insofar as the first 
power of e is concerned. 


Ai(e) 


=£c r 

* j 0 


cos 


2 0 d 0 = 


or 


aio = 1 


The second step is to obtain the first term in the power-series develop- 
ment of Aq (e). Thus, 


Az(e) 



cos 20 cos(0 - e sin 0)d0 
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This expression for A 2 (e) is easily evaluated with the aid of the 
well-known Jacobi expansions in series of Bessel functions; namely. 


00 

cos(e sin 0) = Jo(e) + 2 J2n( e ) cos 2n0 

n=l 


sin(e sin 0) = 2 J 2 n+1 (s)sin( 2 n + l)0 
n=0 


and the power-series expression for Bessel functions of the first kind 


_sa_ 


-2k+n 


j„(«> - 2l (-u k 

k ‘° 2 a+n k:(k + nj: 


Insofar as the second power of e is concerned, the result is 


or 


Similarly, 


or 


a 2( £ ) = “ | e2 


a20 = - 


1 

2 


p r* 

Ao(e) = — e cos(0 - e sin 0)d0 

Jo 


| Ao(e) = i e 2 



and 
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In this manner, the power-series developments for Aq(g) can he system- 
atically constructed. Table I lists values of the coefficients a nm 
sufficient in number to determine the solution, equations ( 3 ), to the 
order . 

Examination of table I shows that the general formula for the set 
of coefficients s^q can be written as 

n n n- l 

a n 0 = -(-l) (n = 1 , 2 , . . • °°) ( 7 ) 

2 n “^- n! 


and that 


^ A2n-i(e) 


n=l 


e 


(8) 


Moreover, from equations (!»•), it can be seen that 


\ Ao(e) +21 Aje) = a 
2 n=l 

Thus, with the aid of equation ( 8 ), it follows that 


77 AqCg) A2n(e) 


( 9 ) 


From considerations of symmetry, the area under the boundary curve 
for one wave length is zeroj that is, 
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Then, by means of equations (ll-). 


r>2rt 

i Ao(e) + An^cos n 0 

CO 

1 - mAn^eJcos m 0 

'0 

_ n=l _ 

hf=1 


d0 = 0 


or 


\ AqCg) = i X! nA n 2 ( £ ) 

2 n=l 


Comparison with equation ( 9 ) shows that 


1 JnV^s) 

2 n=l n=l 


( 10 ) 


Linearized Case 


In the linearized case in which terms involving only the first 
power of e are retained, equations ( 3 ) become 

x = 0 - ee'^sin 0 (lla) 

and 

y = -ijr + ee”^cos 0 (lib) 

It is interesting to note that equation (lla) is similar to one that 
Bessel discussed in connection with Kepler's eccentric-anomaly problem 
(ref. 9)* Corresponding to Bessel's solution. 


00 

0 = x + 2 — J n (nee"^)sin nx (12) 

n=l n 

The process by which this result is obtained gives no information as to 
the conditions under which the expansion on the right-hand side is 
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possible. Therefore, until the interval of convergence has been deter- 
mined, it is not legit im ate to discard terms involving powers of e 
higher than the first. However, the series in equation (12) is a 

Kapteyn series which converges rapidly when ee”^ < 1 and is convergent 

even for ee" 1 ^ = 1. Thus, for the linearized case where e « 1, only 
terms involving the first power of e need be retained. The reversion 
of equations (ll) then yields 


0 = x + ee^sin x 
ijr = y - ee - ycos x 



(13) 


which are in agreement with the results obtained when the linearized 
case is treated in the plane of flow. This result is not in contradic- 
tion with the statements made with regard to equation (l), because for 
the case of vanishingly sma l l, e, the boundary condition of zero normal 
velocity is satisfied along the real axis y = 0 . 


Numerical Example of a Wavy Wall of Finite Amplitude 

In choosing a numerical, example, it is kept in mind that e should 
be less than one but yet should correspond to a wavy wall of moderate 
amplitude. Thus, the value e = O. 7 O corresponds to a wavy wall whose 
amplitude is about one-ninth of the wave length. The conversion of 
results obtained in the 0 i}r-plane to results in the physical flow plane 
is easily made by means of equations (3). Indeed, these equations are 
particularly well-suited for the construction of streamlines, which is 
usually a long and tedious computation. Equation (3a) shows that 0 
is a continuously increasing function of x, such that the effect of 
increasing x by 2rt is to increase 0 by 2 Jt. Moreover, 

A n (e)e -n ^ r sin n 0 is an odd periodic function of x and, hence, 

n=l 

calculations need be made only for the interval 0 5; x ^ n. 

The following formulas are to be utilized for the calculation of 
fluid speed anfl pressure coefficient at the surface of the wavy wall. 
From equation (2), 


dw 

dz 


= u - iv = 


_S2_ 


<T CJ - 


1 - nA n (e)e -n ^cos n0 - i nAj^eJe n ^sin n0 
n=l n=l 
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where u = ^ and v = ^ are the components of the fluid velocity 
ox dy 

in the direction of the x- and y-axes, respectively. 

The magnitude of the fluid velocity at the solid boundary i|r = 0 
is given by 


q 2 = 7= 


1 - 2 nA n (e)cos n0 
n=l 


1 


nAn(e)sin n 0 


n=l 


( 14 ) 


In particular, the maxi miim and ml nimmn speeds occur at 0 = 0 and 
0 = it, respectively; thus. 


1 

( 3max ~5T~ 

1 “ 21 
n=l 


and 


Qmin — 


1 

1 - 21 (~l) n nA n (e) 
n=l 


( 15 a) 


( 15 b) 


The pressure coefficient Cp, obtained by means of Bernoulli's 
theorem, is given by 

Cp = 1 - q 2 ( 16 ) 

The coefficients An(e) have been developed as power series in e 
in the form given by equations (6). Actually, there is no basic reason 
why these coefficients should be expressed in this fashion. In fact, 
this power-series representation of the coefficients An(e), although 
suitable enough for wavy walls of small amplitude (e « l), may not be 
suitable for walls of even moderate amplitude (e < l). Moreover, it 
would appear that for walls of small amplitudes the linearized results 
given by equations (ll) or (13) suffice and that the labor expended in 
obtaining the coefficients is hardly worthwhile . Therefore, in 

order to obtain results appropriate for values of the amplitude e of 
the order of unity, for example, some system of appro xim ation other 
than power-series development in e must be adopted. The most obvious 
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one is that first introduced by Stokes in connection with his work on 
the theory of oscillatory waves (p. 317 of ref. 10 ). Stokes' method 
of approximation is to proceed according to powers of the first coeffi- 
cient Ai(e). That this approach is an independent one can be seen 
from the recursion formula (eq.. ( 5 )) since e appearing therein can be 

OQ 

replaced by the sum Ap^.i (e) (see eq.. (8)). However, the coeffi- 

n=l 

cients already listed in table I can be utilized to obtain the 

expressions for AgCe), A^e), ... as series in ascending powers 

of Ai(e). They are as follows: 


Aa - - i i + A Al 6 - Al 8 ♦ 


A^ = 3*3 _ 53 A 5 _ _Z_ A -,7 + ^^2 A ,9 + 

3 8 Al 384 1 512 Al + 45 x 2-3 Al 


A 1.4. 373 . 6 . 1 . 8 . 

A 4 = - - A! + =■ Al + Al + 

■> 15 X 2 1 9 X 2^ 


= M A, 5 - A,9 + . . . 


Ar = — - Ap 


384 45 X 2 ? 


3 X 2 


,12 


Al" 


> 


(17) 


A6 = .gL Al ^ + J^ l- Al 8 + 

So 315 x 2 10 


_ 16807 . A 7 . Ai 9 + 

' 45 x 2 10 35 X 2 1 5 


Ag . - M Al 8 + . . . 

315 


A 9 


= 531441 .9 . 
35 x 2 1 5 
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From equation (8), It follows tliat 


*1 + 1 A i 3 + Jg A i 5 


+ A 7 + ggS: A 9 

9216 X 45 X 2 14 


+ 


and from equation (9 ) , it 


follows that 


6 (18) 



+ t A ± 


+ Al 6 + H9_ A 8 + 

38I4- 1 3072 1 


(19) 


Note that with the exception of the first coefficient, where the values 
are the same, the coefficients of the power series in Ap shown in 
equations (17) are consistently and substantially less than the corre- 
sponding coefficients of the power series in e. Moreover, for a given 
value of e, the corresponding value of Ap is less. Thus, for e = 0. 
the value of Aj_, obtained from equation (l8), is 0.60. The values of 
the other coefficients calculated by means of equations (17) and (19) 
are as follows: 


A2 = -O.169 

Aij = 0.021 

§*>" °- 220 


Aj = O.O7O 

A£ = -0.010 


% = -0.037 

A j = 0.006 


These values of A n when introduced into equations (4) reproduce 
well the boundary curve y = 0.70 cos xj other streamlines (with > 0) 
will be even more accurately computed because of the presence of the 
factor e“ n 4 r in equations (3 ) . Figure 1 shows the wavy wall for 
e = O.7O and the streamlines which were calculated by means of equa- 


tions (3) for values corresponding to ijr = p, 1, and 2. However, 

when these values of An are inserted into equations (l4) and (16) in 
order to obtain the fluid speed and the pressure coefficient at the 
surface, the resulting points are scattered in a manner which indicates 
a slow convergence of the derived series which appear in these equations 


It is clear that, for wavy walls of even moderate amplitude, 
further analytical development of the coefficients An(e) is of little 
avail because on the one hand a large number of terms of the series 
would have to be taken, and because on the other hand the labor of the 
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approximation increases inordinately with the order of the terms. 
Therefore, some independent method must be found which avoids the 
necessity of determining analytically an infinite number of coefficients. 
This statement may very well apply to flow problems in general when the 
exact solution is to be approximated in some manner by an iterative 
procedure. For the present problem, such a method is described in the 
following section. 


Integral-Equation Approach 
From equation (4b), 


^(e) = — e cos x cos n0 d0 
* Jn 


(n = 0,* 1, 2, ... . »,) 


The evaluation of the infinite number of coefficients A n (e) can be 

made to depend upon a single equation which is obtained by eliminating 
AjjXe) from equation (4a). Thus, 


x = 0 - — e sin n0 / cos x 1 cos n0' d0' 
* n=l J 0 


( 20 ) 


where x ' denotes x as a function of 0 1 . The order of the two 
processes summation and integration cannot immediately be transposed, 
because the series would not be convergent. However, equation (20) 
may be regarded as the limit of the expression which is obtained by 
introducing r n into the general term of the series, where r is a 
quantity less than 1 and in the limit is equal to 1. Thus, summing 
the series 


r^in n0 cos n0' 


leads to 

r sin(0' +0) r sin(0' - 0) 

1 - 2r cos(0' + 0) + r^ 1 - 2r cos(0' - 0) + r^ 




cos x ' d0 1 
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In the limi t 


y 


x = 0 - 



jeot |(0' + 0) - cot |(0' 



x 1 d0* 


and, after integrating hy parts. 


0 - x 



cos x 1 d0' 
cos 0' - cos 0 


( 21 ) 


This is a nonlinear integral equation closely related to the one derived 
hy Theodorsen and Garrick as the basis of their arbitrary-airfoil theory. 
Tn fact, the procedure of solving it is precisely the one devised by 
Theodorsen and Garrick and described in detail in reference 1. Thus, a 
solution is sought, not by development in series according to some 
parameter, but by the following procedure of successive approximations: 

Let 0^ be an approximation of the velocity potential found in 
any manner whatsoever. Inserting this . approximation into the right- 
hand side of equation (21) results in an improved approxi m ation 0k+l* 

In the same manner, further approximations 0^+2* 0k+3> * * * are 
obtained until the iteration ceases, that is, until a pair of succes- 
sive approximations is no longer distinguishable within the limits of 
the desired accuracy. Thus, equation (2l) is written as follows: 


0k+l - x = f sin 




COS X 


* *0k' 


cos 0jj- 1 - cos 0^ 


( 22 ) 


Following the idea of Stokes, 0q is taken equal to x + Aq(e)sin x. 
Than to the second order of approximation. 


01 - x = — sin 

r± it 



cos x 1 d0Q ' 
cos 0o ' - cos 0Q 


( 23 ) 
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Now, according to reference 9 (congiare with eqs. (lla) and (12)), it 
follows that x as a function of 0q is given by 


x = 00 + ^ |(-l) n J n (nA 1 )sin n0 Q 
n=l 


(A X 1) 


and 


cos x = i A-i - 
2 L 


(-l) n 


|j n _l(nAi) - J n+1 (nAi)Jcos n0 o 


Hence, with the aid of the well-known def ini te integral 


r cos n0* d0' _ sin n0 

^0 cos 0' - cos 0 s -*- n 0 


equation (23) yields 


co y \H 

01 - X = -e 53 [jn-l(nAl) - Jn+l^l)] 8 ^ Q0O (2*0 


or 


00 / \H p— 

01 = x - e - ~ n — [Jn-lO^l) - J n+ i(nAiYJ sin n(x + Aq sin x) 


From this point it is not feasible to continue this process in a com- 
pletely analytical manner. For a given value of Ai % 1, however, it 
is a straightforward computation to tabulate corresponding values of 
x and in the range 0 to jt and, by means of a Fourier analysis, 

to express cos x as a coBine series in 0q. (See eq. (25) for means 
of calculating the value of Aq corresponding -to a given value of e.) 
Then, from equation (22), 02 - x is obtained as a sine series in 0q. 
In constrast to the analytical development of the coefficients as power 
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series in e , this method for solving the integral equation (2l) can he 
continued ind efinitely without increase in the amount of computational 
labor for each step. It is clear that, in the limit k — this pro- 
cess yields equation (4a) hut, in practice, it is terminated when the 
difference between two successive approximations is deemed negligible. 

For values of Ap(e) greater than unity, the appro xima tion 
00 = x + Ap(e)sin x is not valid because x then becomes a multiple 
root of this equation in the range 0 ^ x % it. The requirement is that 
0 be a continuously increasing function of x within this range at 
every stage of the iteration procedure. For Ap > 1, a plausible 
beginning would be 0q = x + sin xj it may be presumed that this initial 
relationship between x and 0 O leads to a convergent process up to a 
certain value of e. The last valid relationship between x and 0 
would then agal n be taken as the beginning of the iteration procedure 
for a still higher range of values of e. This process of continuation 
may presumably be carried on indefinitely. 

Now, for some range of values 0 < Ap < 1 (but possibly for values 
of e > l), equation (24) may be considered accurate enough to terminate 
the process of iteration. That is, approximately. 


x = 0 + e ^ jj n _ 1 (nA 1 ) - J n+1 (nAi)Jsin n0 

Comparison of this equation with equation (4a) shows that 
An ( € ) = _ e iz^jj^nAp) - J n+ p(nApj] 


(25) 


Note that for n = 1, 


Ap(e) = e[j 0 (Ap) - J 2 (Ap)] 


which provides a transcendental equation for the determination of Ap 
for a given value of e. With the use of the recurrence relation 


Jn-i(nAi) - J rt+-p(nAp) 


2 dJnCnAi) 
n 
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it follows that 


•An( e ) 


(-l) n dJn(nAi) 
n 2 ‘ dAi 


Then 



d y- J2n-1 [fo - 

dAj_ n=l (2n - l)^ 


However, 


00 



Jgn-igai - DaJ 

(2n - l) 2 



and hence (see eq^. (8)) 


n=l 


Agn-l^ 6 ) 


e 


Again (see eq[. ( 9 )) 


\ Ao(e) ^ A2n( e ) 

^ n=l 


or 


\ Ao(s) = \ e 

d 2 dAi n=1 


Jpn(2nAi ) 
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But 


J 2n( 2nA l^ _ 1 a ^2 


n=l n c 


and hence. 


| Ao(e) = \ eA X 


Then to the second order of approximation, equations (3) may he •written 
as 


00 n 

x = 0 + e 21 - ^ - } --[j n _p(nAp) - J n+ p(nAp)] e _n tsin n0 

n=l *— - J 


00 

y = 4r + | eAp - e 2Z jj n -l(nAp) - J n+ p(nAp)j e^cos n0 (Ap < l) 

(26) 

OO 

It is of interest to examine the sum S = 2 nAj^eJsin n0 which 

n=l 

appears in equation ( ih- ) for q^. To the second order of approximation. 


S = -e ^>2 (-l^l^n-l^l) " Jn+lC^All sin n 0 
n=l L 


Now, hy defining a new variable p hy the equation 


(27) 


0 = p + Ap sin p 


then 


00 / >11 

p = 0 + 2 — J n (nAp)sin n0 


n=l 
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and 


cos 


00 , » n 

(i = | Ai - ^ — [Jn-l^l) ~ Jn+iCnAiTjcos n0 


Hence, 


sin 


H- ^ ^ (-l) n jjn-l(nAi) - Jn+i(nAq7J sin n0 


or, since 


dp. _ 1 

d0 1 + Aq cos p 


S = 


e sin p 
1 + Aq cos p 


(Ai< lj 0< p < it) 


The graph of S as a function of 0 for the case e = O.7O (and 
Ai = O.606 from eq. (25)) is shown as the solid curve in figure 2. 

The circles correspond to the evaluated series, equation (27), with 
n = 1 to 24. The values of A n according to the second order of 
approximation are 


Aq = 0.606 

A 2 = -0.164 

A3 = O.067 

A4 = -0.033 

Ar> = 0.017 

A£ = -0.010 

A j = 0.006 

\ Aq = 0.212 



These values agree remarkably well with those obtained by the method of 
Stokes. Note, however, that only one step was needed in the solution 
of integral equation (2l) to achieve the accuracy of nine increasi ng ly 
laborious steps according to powers of the first coefficient Aq(e). 

The plus signs in figure 2 correspond to points computed with the Stokes 
values of A^. The scattering of these points, particularly in the 
critical portion of the curve, is due to the slow convergence of the 
series for S with only the first few terms available and presages 
inaccurate determination of the derived quantities, fluid speed and 
pressure coefficient at the surface of the solid boundary. Here again 
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the superiority of the integral- equation approach is demonstrated in 
that the second-order approximation provides an analytic form for the 
general coeffficient An(e) (see eq. (25))- Table II lists the results 
of the calculation of fluid speed and pressure coefficient at the sur- 
face of the wavy wall to the second order of approximation, and figure 3 
shows the correspo nding graphs for an interval of one wave length. 

Tn conclusion, the analysis of the present paper indicates the 
general unde sir ability of determining the velocity potential or the 
stream function in the form of a series whose coefficients are ana- 
lytically developed in a power series of some parameter. Rather, the 
flow problem should be set up in the form of an integral equation whose 
solution can be obtained by a method of successive approximations. If 
this course is not possible, the parameter for the analytical develop- 
ment of the coefficients should be chosen so as to insure rapid con- 
vergence of the series. As the example of the present paper shows, 
this parameter may not be the most obvious one. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 25, 1953* 
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TABES I.- COEFFICIENTS 


n 

a nm 

m = 0 

m = 1 

m = 2 

m = 3 

warn 

0 


1 

1 

11 

143 


2 

8 

192 

46o8 

1 

1 

3 

15. 

1471 

82649 

8 

64 

9 X 2 10 

45 X 2 1 ^ 

2 

HlCVJ 

1 

ll 

24 

_ 53 _ 
128 

16649 
45 x 2 10 


3 

3 

215 

683 

262801 


8 

384 

1024 

45 x 2 15 



1 

1333 

11873 



*T 

" 3 

15 x 2 7 

45 X 2 8 



5 

125 

384 

40187 
45 x 2 10 

190849 
15 X 2 15 



6 

27 

356983 





80 

315 x 2 10 




7 

7 5 

4870981 




45 x 2 10 

105 X 2 1 5 




8 

128 

315 





9 

9 6 





35 x 2 1 5 
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TABLE II— CAILJULATED VALUES OF FLUID SPEED AND PRESSURE COEFFICIENT 


0 

B 

S 

24 

y An sin n0 
n=I 

24 

~ An cos n0 
n=°l 

24 _ 

nAn sin n0 

n=I 

24 

y nAn cob n0 

X 

<1 

Op 

0 

0 

0 

0 

0.4879 

0 

0 . 4 oi 

0 

I.67O 

-1.789 

Jt 

4 

.496 

.217 

.3063 

.4033 

.218 

•367 

.479 

1.494 

- 1.232 

it 

2 

1.046 

.465 

.5523 

.1385 

.465 

.238 

1.019 

1.120 

-.254 


1.490 

.665 

.6300 

-.1559 

.665 

.034 

1.464 

.852 

• 273 


1.760 

.776 

.6162 

-.3445 

.776 

-.153 

1.740 

.720 

.482 

5 * 

6 ' 

0 . 864 jt 

2.093 

2.222 

.870 

.880 

,•5177 

-.5610 

.869 

-.465 

2.100 

.587 

.656 








rut 

12 

2-554 

.795 

.3515 

-.7867 

•795 

-1.019 

2.529 

.461 

.788 

to 

48 

2.977 

.285 

.1004 

-.9026' 

.285 

-1.501 

2.976 

• 397 

.842 

jt 

it 

0 

0 

-.9121 

0 

-1.550 

rt 

•392 

.846 
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Figure 2 
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Figure 3 •- Fluid speed and pressure coefficient at surface of wavy wall, 

y = 0.70 cos x. 
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